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Abstract. Nowadays, the analysis of the X-ray spectra of magnetically powered neutron stars 
or magnetars is one of the most valuable tools to gain insight into the physical processes 
occurring in their interiors and magnetospheres. In particular, the magnetospheric plasma 
leaves a strong imprint on the observed X-ray spectrum by means of Compton up-scattering 
of the thermal radiation coming from the star surface. Motivated by the increased quality 
of the observational data, much theoretical work has been devoted to develop Monte Carlo 
(MC) codes that incorporate the effects of resonant Compton scattering (RCS) in the modeling 
of radiative transfer of photons through the magnetosphere. The two key ingredients in this 
simulations are the kinetic plasma properties and the magnetic field (MF) configuration. The 
MF geometry is expected to be complex, but up to now only mathematically simple solutions 
(self-similar solutions) have been employed. In this work, we discuss the effects of new, more 
realistic, MF geometries on synthetic spectra. We use new force-free solutions [14] in a previously 
developed MC code to assess the influence of MF geometry on the emerging spectra. Our 
main result is that the shape of the final spectrum is mostly sensitive to uncertain parameters 
of the magnetospheric plasma, but the MF geometry plays an important role on the angle- 
dependence of the spectra. 



1. Introduction 

Neutron stars (NSs) can be powered by different sources of energy: rotation, accretion, residual 
heat and magnetic field. Up to now, more than twenty X-ray pulsating sources are identified 
as isolated, slowly rotating (P ~ 2 — 10 s) neutron stars, whose large X-ray luminosity 
Lx ~ 10^^ — 10^^ erg/s cannot usually be explained in terms of rotational energy, at variance 
with rotation-powered pulsars (RPPs). Historically, they are classified as Anomalous X-ray 
Pulsars (AXPs) or Soft Gamma Repeaters (SGRs), and they are recognized as magnetars. In 
these objects, MF decay and evolution is the dominant energy source and crustal deformations 
produced by the huge internal magnetic stresses are the likely cause of bursts and episodic giant 
flares. In most cases, the spin-down inferred external magnetic field turns out to be super-strong. 
Recently, increasing evidence gathered in favor of a unified scenario for isolated neutron stars. 



their different observed manifestations being mainly due to a different topology and strength of 
the internal MF, and to evolution 

The typical soft X-ray (< 10 keV) spectrum of magnetars is given by a thermal component 
(kT ~ 0.3 — 0.5 keV), plus a hard tail phenomenologically described by a power law with a 
photon index F ~ 3 — 4 [7]. This hard tail is believed to be produced by resonant compton 
(up) scattering (RCS) of photons emitted by the star surface onto mildly relativistic particles 
flowing in magnetosphere. In order to take into account for this effect, several codes have been 
developed in recent years [HEHSJE]. The two basic inputs are the external MF and a prescription 
for the particle space and velocity distribution. With these ingredients, and with a given energy 
distribution for the incoming surface photons (the seed spectrum), the MC simulations of RCS 
provide the emerging, reprocessed spectrum, which must be compared with observational data. 

Magnetars also exhibit a high-energy (~ 20-200 keV) power-law tail, the origin of which is 
still under debate. One of the possible mechanisms is again RCS, whose efficiency at high energy 
is very sensitive to MF and plasma velocity. In this paper we discuss the effect of the geometry 
of the MF, which is obtained from new calculations of magnetostatic, force-free equilibria. 

2. Force-free tw^isted magnetosphere 

The shape of the internal MF of NSs is largely unknown, but it is thought to be complex, with 
toroidal and poloidal components of similar strength, and with significant deviations from a 
dipolar geometry. In the crust, a large toroidal field is a necessary ingredient to explain the 
observed thermal and bursting properties of magnetars, even in the case in which they show 
a low value of the inferred external dipolar field [13] . The complex internal topology must be 
smoothly matched to an external solution, so that the internal MF geometry at the star surface 
has a strong effect on the external MF. In this framework, the RPP phenomenolgy is expected 
to be better explained by the presence of a NS with a simple, nearly potential, dipolar field, 
while magnetars activity is compatible with the presence of a complex external field, as a result 
of the transfer of magnetic helicity from the internal to the external field. A non-potential MF 
needs supporting currents, and charges can be either extracted from the surface or produced by 
one-photon pair creation in the strong MF [1]. 

In the inner magnetospheric region MF lines corotate with the star, with footprints connected 
with the crustal field, which evolves on long timescales (10^ — 10^ years). Furthermore, resistive 
processes in the magnetosphere act on a typical timescale of years [2], much longer than the 
typical response of the tenuous plasma, whose Alfven velocity is near to the velocity of light. 
For these reasons, a reasonable approximation is to build stationary, force- free magnetospheric 
solutions, assuming that the footprints are anchored in the crust and neglecting resistivity. Then, 
the ideal MHD force-free equation is a simple balance between electromagnetic forces, provided 
that inertial, pressure, and gravity terms are negligible 

PeE + -JxB = 0, (1) 
c 

where we have introduced the rotationally induced electric field E = —Vrot x B, the related 
charge density and the current density J = cV x B/4tt. 

The electric field gives corrections of order {pQ./cf' . Neglecting such contributions is safe if 
we consider the region closest to a magnetar surface (r < lOOr^), and remembering that the 
light cylinder for slow rotators lies typically at a few thousands stellar radii [ri = c/il and 
magnetar periods are in the 1-10 s range). Therefore, condition ([1]) reduces simply to J || B: 
the particles can flow only along MF lines. Even if the pulsar equation simplifies in this limit, 
it still contains an arbitrary function, which basically describes the toroidal field as a function 
of the poloidal one (see for the mathematical formulation). In the literature, only the 
simplest (semi-)analytical solutions have been derived and later employed in applications. In 
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Figure 1. Numerical model 1 (left) and model 2 (right). First panel: \B\ (grey logarithmic 
scale) with scattering surfaces for photons of 1, 3 and 5 keV (here /3 = 0.5, n = 0.5; see text); 
second panel: current intensity \ J\/c (grey linear scale), in units of 10^^ G/r^. 



particular, in the context of magnetars, the self-similar twisted dipole family of solutions is the 
most popular choice. Proposed for solar corona studies [6], and later extended to the magnetar 
scenario [,1 2J , in self-similar models the radial dependence of the MF components is described by 
a simple power law, Bi oc r~^P~^'^\ These solutions are semi-analytical and described by only one 
parameter, the radial index p (see [10] for globally twisted multipolar fields). However, globally 
twisted solutions necessary evolve in time and magnetic configurations in which the twist affects 
only a limited bundle of field lines are likely required by observations [2j . 

As an alternative to self-similar models, in [T3] we presented a numerical method to build 
general solutions. The force- free solutions are obtained by iteration starting from an arbitrary, 
non-force-free poloidal plus toroidal field. By the introduction of a fictitious velocity field 
proportional to the Lorentz force, the code is designed to dissipate only the part of the current 
not parallel to the MF. The boundary conditions are an arbitrary, fixed radial MF at the surface, 
i?r(r^, 0), and the requirement of a continuous matching with a potential dipole at large distances 
''"out > lOOr^,. This method is able to produce general configurations, whose features may be 
radically different from the self-similar solutions. In this work, we try to reproduce some j- 
bundle like structure, aiming at providing input for MC simulations [9]. In particular, we focus 
on two models, whose MF and current distribution are shown in Fig. [TJ The first one is obtained 
fixing a dipolar form for radial MF, while the second one has a concentration of field lines near 
the south pole. To get these profiles, we have chosen the (j) component of the vector potential 
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from which we recover 5^ = (V x A^t/)),- (-Bo is the value of Bj. at north pole). Both models 
have high helicity, but it is not homogeneously distributed like in self-similar models. Model 1 
is characterized by an equatorial asymmetry, with current and twist concentrated in a closed 
bundle near the equatorial region and a j— bundle near the southern semi-axis. Model 2 is more 
extreme: its helicity is ~ 25 times larger than for model 1. Both are supposed to be more 
realistic than the self-similar models, since currents are more concentrated near the axis. 



3. Resonant Compton scattering 

The MC code we employ is presented in [HI [9] and we refer to these works for further details. 
Here we just summarize some basic concepts. It is assumed that currents flow in a non-potential 
magnetosphere and that the electron density is large enough for the medium to be optically thick 
to resonant scattering: photons produced by the cooling star surface can be then up-scattered 
by RCS, populating the hard tail of the spectrum. In the stellar frame, the resonant energy 
corresponds to the Doppler-shifted gyration frequency of the electron 

u;DiB,P) = —l—— (4) 
7(1 - ruec 



where B = \B\ is the local intensity of MF, c/3 the electron velocity, 7 = — the Lorentz 

factor, and n the cosine of the angle between the photon direction and the electron momentum. 
The non-resonant contributions are negligible, since the resonant cross-section is up to five orders 
of magnitude larger than the Thomson cross-section. Figs. [T]show the scattering surfaces in the 
two models, for photons of 1, 3 and 5 keV, assuming /3 = 0.5 and fi = 0.5. The surfaces are far 
from being spherically symmetric. In model 2, which has the strongest helicity, the scattering 
surfaces of soft X-ray photons lie tens of stellar radii away from the surface, especially near the 
axis. 

The scattering optical depth for a given seed photon energy depends on the local particle 
density Tig, the velocity of the scatterer and the local value of the MF. Since photons can be 
pushed to energies of some hundreds keVs, non-conservative scattering and the fully relativistic 
QED cross sections are used [9]. The velocity distribution is assumed to be a one-dimensional 
(along the field line), relativistic Maxwellian distribution with a plasma temperature Tg and a 
bulk velocity Phuik, as described in [8j. The strong simplification in this approach is that the 
velocity distribution does not depend on position, which reduces the microphysical inputs to 
two parameters: the temperature of the plasma and Pbuik- Only recently, a more consistent 
treatment of the charges space and velocity distribution was attempted ^ . These new results 
show that the velocity distribution is non-trivial and strongly depends on position, which opens 
a new range of possibilities, so far unexplored. 

Given a seed spectrum for the photons emitted from the surface (assumed as a blackbody 
with temperature T^), the MC code follows the photon propagation through the scattering 
magnetosphere. When a photon enters in a parameter region where no more resonant scatterings 
are possible, the photon escapes to infinity, where its energy and direction are stored. The sky at 
infinity is divided into patches, so that viewing effects can be accounted for (e.g. if line-of-sight 
is along the angles Og, (pg, only photons collected in the patch which contains that pair of angles 
are analyzed). The angle- averaged spectrum is obtained simply by averaging over all patches. 



4. Results and discussion 

In previous works [15], self-similar magnetospheric models were used to obtain synthetic spectra 
that were satisfactorily fitted to X-ray observations, giving typical values of /3huik ^ [0-1, 0.7] and 
A(j) S [0.4, 2.0]. In this paper we explore the dependence of the spectrum on the magnetosphere 
model by comparing results obtained with self-similar models and with our new numerical 
solutions for force-free magnetospheres. We present results comparing models with a magnetic 
field intensity at the pole of Bp = 10^^ G. The seed spectrum is assumed to be a kT^, = 0.4 keV 
Planckian isotropic distribution for both ordinary and extraordinary photons. In Figs. [2] and 
[3] we plot the angle-averaged spectra, i.e., the distribution of all photons escaped to infinity, 
independently on the final direction. Fig. [2] shows the effects produced by changes in the 
bulk velocity, keeping fixed the MF configuration (model 1), and the electron temperature 
(kTe = 20 keV). The main fact is the the spectrum in the region > 10 keV becomes harder as 
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Figure 2. Angle-averaged spectra for kTg = Figure 3. Angle- averaged spectra for kTe 

20 keV, model 1, changing /Shuik- 20 keV, /3buik = 0.9, changing MF model. 



I^buik is increased, while it is depleted of photons of energies in the 1 — 10 keV range. A similar, 
but less pronounced effect, is obtained by increasing the electron temperature [9]. 

In Fig. [3] we compare spectra obtained with fixed values of kT^ = 20 keV and f^huik = 0-9, 
varying only the MF topology. The three lines correspond to model 1 (solid), model 2 (dashes) 
and a self-similar model (dash-dotted line) with A(/> = 1.36, which approximately has the same 
total helicity of model 1. We have chosen a high value of f3i,uik to show a case where the effects are 
larger, but our conclusions do not change qualitatively for lower values of (3buik- The high-energy 
cutoff produced by electron recoil is also evident. We see that the major differences arising from 
changes in MF topology appear in the hard tail of the spectra, which is clearly harder for the 
self-similar model. However, the thermal part of the spectrum is more depleted, especially for 
model 2. 

However, in the angle-averaged spectrum, the differences due to the MF topology are partially 
smeared out, and the comparison between synthetic spectra as seen from different viewing 
angles (first three panics in Fig. U]) is more interesting. The comptonization degree reflects the 
inhomogeneous particle density distribution and therefore the particular current distribution 
produces important differences between different viewing angles. Since neutron stars rotate, 
the study of pulse profiles and phase-resolved spectra can trace the geometric features of the 
scattering region (tens of stellar radii). In the lower right panel we show the light curves for the 
three MF models, in bands 0.5-10 keV and 20-200 keV, for an oblique rotator with a certain 
hue of sight (^ = x = 30° in [5]). 

For self-similar models (lower left panel), the thermal part of the spectrum shows a smooth 
dependence with the viewing angle, which in turn translates into relatively regular light curves. 
In models 1 and 2, differences in the spectra induced by viewing angles are more pronounced. 
In these cases, the spectrum is much more irregular and asymmetries are larger. In particular, 
model 1 (upper left) shows a softer spectrum when seen from northern latitudes, with important 
spectral differences. The pulsed fraction of model 1 is very high in the hard X-ray band, while 
it is comparable with the self-similar model in the soft range. Model 2 (upper right) has a more 
symmetric distribution of currents, and the comptonization degree at different angles depends 
on the energy band in a non trivial way. This results in large pulsed fraction for both energy 
ranges and in notable differences between their pulse profiles. 

The high energy tail, seen at different colatitudes, can vary by one order of magnitude or 
more. Investigating the geometry can help to recognize the different components seen in the hard 
tails via pulse phase spectroscopy. On the other hand, the macrophysical approach used in this 
work to obtain the MF topology should be accompanied by the corresponding microphysical 




Figure 4. Synthetic spectra for model 1 (upper left panel), model 2 (upper right), self-similar 

model (lower left) with the parameters indicated in the figures. Legends indicate the four viewing 
angles here considered: aligned with north or south pole, just above and just below the equator. 
The seed blackbody is shown for comparison (solid line). Fourth panel shows the pulse profiles 
obtained for model 1 (solid line), model 2 (dotted), self-similar (dashed), in the energy ranges 
0.5-10 keV (thin lines) and 20-200 keV (thick lines), with ^ = x = 30° (see text). 



description to determine the velocity distribution of particles, here simply taken as a model 
parameter. A fully consistent description of both MF geometry and the spatial and velocity 
distribution of particles is needed to advance our interpretation of magnetar spectra. 
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